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Abstract 

In  this  manuscript,  embedded  atom  potentials  are  used  to  determine  the  various  core  structures  of  glide  screw-character  dislocations 
intersecting  free  surfaces  at  right  and  inclined  angles  in  face-centered  cubic  Ni  and  Cu.  It  is  shown  that  the  negative  constriction  forms  at 
free  surfaces  under  certain  conditions.  The  role  of  various  factors  affecting  the  formation  of  negative  constriction  at  the  free  surface, 
negative  constriction  energy,  ledge  annihilation  and  screw  dislocation  rotation  due  to  Lothe’s  forces  are  discussed.  The  activation  energy 
for  surface  cross-slip  nucleation  when  the  screw  dislocation  intersects  the  free  surfaces  at  right  angles,  in  the  absence  of  ledge  annihilation 
forces,  is  shown  to  be  0.05  eV  in  Ni  and  0.09  eV  in  Cu.  The  activation-energy  values  obtained  via  the  nudged  elastic  band  method  are 
significantly  lower  than  the  activation  energy  for  cross-slip  at  attractive  forest  dislocation  intersections  in  these  materials.  The  present 
results  are  expected  to  be  useful  in  understanding  the  mechanical  behavior  of  micron-sized  pillars  under  pure  tension  and  compression. 
©  2013  Acta  Materialia  Inc.  Published  by  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

There  is  increasing  recognition  of  the  need  to  incorpo¬ 
rate  physics-based  models  of  deformation  in  the  design  of 
structural  components.  While  models  for  predicting  yield 
strength  and  creep  behavior  are  beginning  to  incorporate 
significant  physics,  models  for  fatigue,  ultimate  strength 
and  ductility  remain  mostly  empirical.  This  is  because 
strain-hardening  and  fatigue  resistance  are  strongly  influ¬ 
enced  by  dislocation  micromechanisms,  especially  cross¬ 
slip,  and  including  physics-based  cross-slip  processes  in 
mesoscale  simulations  has  been  difficult.  The  early  work 
of  Escaig  remains  the  most  widely  cited  and  used  model 
for  cross-slip  wherein  the  screw  dislocation  locally  con¬ 
stricts  and  immediately  redissociates  on  the  cross-slip  plane 
[1];  however,  this  model  poses  several  difficulties  with 
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respect  to  quantitative  simulations.  The  model  is  highly 
sensitive  to  choice  of  parameters  that  have  thus  far  been 
difficult  to  quantify.  For  example,  the  constriction  width 
required  for  cross-slip  is  unknown  but  significantly  influ¬ 
ences  the  energetics  of  the  cross-slip  process  [2].  This  diffi¬ 
culty  has  led  to  an  ad  hoc  postulate  that  obstacles  always 
exist  in  materials  and  that  they  enable  sufficient  dislocation 
core  constriction  under  normal  stresses,  thereby  ensuring 
cross-slip  [3].  Aside  from  being  unsatisfactory,  this  forces 
cross-slip  models  (particularly  for  single  crystals)  to  make 
arbitrary  assumptions  about  dislocation  obstacle  spacings. 
Advances  in  atomistic  simulations  offer  the  possibility  of 
gaining  insights  into  the  cross-slip  process  that  may  serve 
to  enable  mesoscale  simulations  to  accurately  capture  the 
atomic-level  physics  of  a  dislocation  process. 

Atomistic  simulation  studies  from  the  late  1990s  were 
limited  to  calculating  constriction  energies  and  energetics 
of  cross-slip  using  a  single  dislocation  in  a  periodic  unit  cell 
[2,4].  These  simulations  helped  to  unravel  some  interesting 
and  useful  aspects  of  cross-slip  which  were  unanticipated. 
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The  most  important  finding  was  that  configurations  exist 
for  which  one  of  the  pair  of  constrictions  has  a  negative 
energy  [2,4].  Thus  spontaneous  nucleation  of  cross-slip  is 
possible  so  long  as  only  one  constriction  is  needed,  as  is 
the  case  for  surface  nucleation.  However,  previous  atomis¬ 
tic  simulations  [4]  of  a  screw-character  dislocation  inter¬ 
secting  free  surfaces  at  right  angles  in  face-centered  cubic 
(fee)  Cu  have  shown  that  the  screw  dislocation  prefers  to 
reside  on  a  {11  l}-type  plane  parallel  to  the  surface  ledge 
on  which  it  was  created,  without  easy  cross-slip  nucleation. 
This  is  in  contrast  to  recent  atomistic  simulation  results  of 
screw-character  dislocations  intersecting  mildly  attractive 
and  repulsive  forest  dislocations  in  the  bulk  in  fee  Ni, 
where  cross-slip  has  been  shown  to  occur  easily  [5-8].  In 
this  work,  the  possibility  of  cross-slip  at  dislocation-sur¬ 
face  intersections  in  fee  Ni  and  Cu  was  examined  using 
large-scale  atomistic  simulations.  A  glide  screw  dislocation 
that  intersects  two  free  surfaces  perpendicularly  as  well  as 
in  an  inclined  fashion  is  modeled  and  it  is  shown  that  sev¬ 
eral  conditions  exist  where  surface  cross-slip  can  nucleate 
easily.  This,  we  believe,  has  the  potential  to  explain  disloca¬ 
tion  multiplication  in  low-density  micron-sized  fee  pillars 
subjected  to  a  uniform  tensile  or  compressive  loading. 

We  have  performed  molecular  statics  simulations  of 
screw  dislocation  intersections  with  free  surfaces  for  the 
[111]  and  [100]  orientations,  with  the  screw  dislocation 
having  a  Burgers  vector  of  1/2  [iTO].  Section  2  describes 
the  simulation  technique  and  the  interatomic  potentials 
used  for  the  simulations.  Section  3  presents  the  core  struc¬ 
tures  obtained  from  atomistic  simulations.  Section  4 
describes  activation  energy  calculations  for  surface  cross¬ 
slip  when  the  screw  dislocation  intersects  the  free  surfaces 
at  right  angles,  in  the  absence  of  ledge  annihilation  forces. 
Finally,  Sections  5  and  6  give  a  discussion  and  summary  of 
the  results,  respectively. 

2.  Simulation  technique 

The  atomistic  simulations  described  here  employed  the 
three-dimensional  (3-D)  parallel  molecular  dynamics  code, 
LAMMPS  [9],  developed  at  Sandia  National  Laboratory. 
The  simulation  cell  is  a  rectangular  parallelepiped,  with 
the  z-axis  oriented  along  the  [111]  or  [100]  directions. 
For  the  z-axis  oriented  along  [111],  the  dimensions  of  the 
simulation  cell  along  the  z-axis  was  ~  1 5  nm,  the  dimension 
along  the  x-axis  was  ~30.0  nm  and  along  the  y-axis 
~15nm.  For  the  z-axis  oriented  along  [100],  the  dimen¬ 
sions  of  the  simulation  cell  along  the  z-axis  was  ^40  nm, 
the  dimension  along  the  x-axis  was  ~14nm  and  along 
the  y  axis  ~14nm.  This  corresponds  to  approximately  1 
million  atoms  in  the  simulation  cell.  A  1/2  [110]  screw  dis¬ 
location  is  inserted  in  the  middle  of  the  simulation  cell 
using  its  anisotropic  elastic  displacement  field  [10].  Unless 
otherwise  specified,  the  plane  of  the  cut  performed  to  gen¬ 
erate  the  screw  dislocation  was  the  (111)  glide  plane.  In  the 
activation  energy  calculations,  the  cube  plane  cut  was  used. 
The  origin  for  the  initial  anisotropic  displacement  field  of 


the  screw  dislocation  was  varied  to  obtain  several  different 
core  structures  for  the  screw  dislocation  interaction  with 
the  free  surfaces.  Fixed  boundary  conditions  were  applied 
along  the  z  and  y  directions  and  energy  minimization  was 
performed  using  the  conjugate  gradient  technique.  Free 
surface  boundary  conditions  were  applied  along  the  x 
direction. 

2.7.  Interatomic  potential 

The  embedded  atom  potential  used  in  the  simulations  is 
the  potential  developed  for  fee  Ni  by  Angelo,  Moody  and 
Baskes  [11]  (hereafter  termed  the  Moody  potential)  based 
on  the  Voter  and  Chen  format  as  well  as  the  embedded 
atom  potential  developed  for  fee  Ni  by  Mishin  [12].  For 
fee  Cu,  an  embedded  atom  potential  developed  by  Mishin 
et  al.  is  used  [13].  Table  1  gives  the  lattice  parameter,  cohe¬ 
sive  energy,  elastic  constants  and  stacking  fault  energy 
given  by  the  Mishin  and  Moody  potentials  for  fee  Ni  as 
well  as  the  Mishin  potential  for  fee  Cu. 

2.2.  Depiction  of  core  structures 

In  some  figures,  to  depict  core  structures  we  use  the 
method  developed  by  Stukowski  and  Albe  [14]  which 
extracts  dislocation  lines  and  their  associated  Burgers  vec¬ 
tors  from  3-D  atomistic  simulations.  It  is  based  on  a  fully 
automated  Burgers  circuit  analysis,  which  locates  disloca¬ 
tion  cores  and  determines  their  Burgers  vector.  The  transi¬ 
tion  from  the  atomistic  system  to  a  discrete  dislocation 
representation  is  achieved  through  a  subsequent  vectoriza- 
tion  step. 

In  order  to  illustrate  the  relaxed  screw  dislocation  geom¬ 
etries,  in  some  figures,  we  take  advantage  of  the  increase  in 
atomic  energy  produced  by  the  strain  field  of  the  partial 
dislocations.  By  plotting  the  atoms  with  assigned  energies 
within  LAMMPs  of  greater  than  —4.42  eV  (Ni)  or 
3.52  eV  (Cu)  (the  energy  of  atoms  in  the  stacking  fault 
region),  the  partial  dislocations  can  be  imaged  easily  even 
in  these  large  simulation  cells  [5-7].  In  order  to  illustrate 
the  cross-slipped-segment  products  of  the  screw  disloca¬ 
tion,  the  positions  are  shown  in  a  [1 1 2  ]  projection  as  well 
as  the  [111]  projection.  In  the  [112]  projection  segments 
spread  on  the  initial  (111)  plane  appear  as  a  single  line 
and  cross-slipped  segments  (i.e.  on  a  (1 1  T)  plane)  appear 

Table  1 


Lattice  parameter,  gq,  cohesive  energy,  Ec,  elastic  constants,  Cu,  C\2  and 
C44,  and  stacking  fault  energy,  7,  given  by  the  Moody  embedded  atom 
potential  for  Ni,  and  the  Mishin  embedded  atom  potentials  for  Ni  and  Cu. 


Property 

Ni  (Mishin) 

Ni  (Moody) 

Cu  (Mishin) 

flo  (A) 

3.52 

3.52 

3.615 

Ec  (eV) 

-4.45 

-4.45 

-3.54 

C„  (xl0"Nm-^) 

2.41 

2.46 

1.70 

Ci2(xlO"Nm-^) 

1.51 

1.47 

1.23 

C44  (xlO"  Nm“^) 

1.21 

1.25 

0.76 

y  (mJ 

134 

89 

44 

2 
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as  a  pair  of  partials  separated  by  a  stacking  fault.  This 
method  of  depicting  the  core  of  the  screw  dislocation  also 
readily  shows  the  sense  of  rotation  of  the  screw  dislocation 
at  the  free  surfaces. 

3.  Screw  dislocation  intersecting  the  free  surfaces  at  right 
angles 

Figs.  1  and  2  give  a  dislocation  line  representation  of  the 
core  structure  obtained  for  the  screw  dislocation-surface 
interaction  when  the  screw  dislocation  intersects  the  free 
surfaces  at  right  angles  using  the  Mishin  potential  for  Ni. 
The  z-axis  is  along  [111],  the  x-  and  y-axes  are  along 
[110]  and  [112],  respectively.  In  Figs.  1  and  2,  a  1/2 
[1 T  0]  screw  dislocation,  passing  approximately  through 
the  origin,  was  initially  inserted  into  the  simulation  cell 
using  its  anisotropic  elasticity  displacement  field,  with  sev¬ 
eral  slightly  different  elastic  centers  for  the  displacement 
field.  The  cut  for  the  initial  elastic  displacement  field  was 
made  on  the  [111]  glide  plane.  Depending  upon  the  elastic 
center  position,  either  the  core  structure  shown  in  Fig.  1  or 
that  shown  in  Fig.  2  was  obtained,  having  the  screw  dislo¬ 
cation  either  fully  residing  on  the  (1 1 1)  glide  plane  (Fig.  1) 
or  residing  partially  on  the  (1 1 1)  glide  plane  and  partially 
on  the  (1 1  T)  cross-slip  plane  (Fig.  2).  The  core  structure 
with  the  screw  dislocation  completely  on  the  (1 1  T)  cross¬ 
slip  plane  never  materialized.  Similar  results  were  obtained 
with  the  Moody  potential  for  Ni  as  well  as  the  Mishin 
potential  for  Cu.  In  Fig.  1,  the  screw  dislocation  is  con¬ 
stricted  at  the  positive  ‘x’  surface,  whereas  it  is  expanded 
at  the  negative  ‘x’  surface.  This  is  due  to  the  effect  of  Escaig 
stress  components  acting  on  the  Shockley  partials  arising 
from  image  stresses  at  the  free  surfaces  and  has  been  pre¬ 
dicted  using  continuum  elasticity  theory  [15-17].  In  addi¬ 
tion,  in  Fig.  2,  the  screw  dislocation  is  bent  inwards 
(toward  the  positive  direction)  at  both  free  surfaces. 
This,  we  believe,  results  in  an  annihilation  of  ledge  surface 
area  at  the  free  surfaces  since  the  screw  dislocation  was 
formed  by  performing  a  (1 1 1)  plane  cut.  The  ledge  annihi¬ 


lation  forces  are  balanced  by  line  tension  forces  at  equilib¬ 
rium  and  this  is  the  reason  why  the  screw  dislocation 
prefers  to  reside  on  the  plane  on  which  the  initial  cut  was 
performed.  To  interrogate  the  screw  dislocation  behavior 
at  free  surfaces  without  the  ledge  annihilation  effect,  several 
simulations  having  slightly  different  initial  elastic  centers 
were  performed,  with  the  screw  dislocation  formed  using 
a  (00 1)  cube  plane  cut  instead  of  a  (1 1 1)  plane  cut,  again 
using  the  Mishin  potential  for  Ni.  In  this  case,  four  differ¬ 
ent  core  structures  were  obtained  for  the  screw  dislocation 
as  shown  in  Fig.  3a-d  in  a  3-D  dislocation  line  representa¬ 
tion.  For  two  of  the  four  different  core  structures,  the  screw 
dislocation  is  partially  cross-slipped  via  the  spontaneous 
formation  of  a  negative  constriction  [2]  and  the  two  core 
structures  are  equivalent  by  symmetry.  The  other  two  core 
structures,  one  residing  fully  on  the  glide  plane  and  one 
fully  residing  on  the  cross-slip  plane,  are  also  crystallo- 
graphically  equivalent  to  each  other.  Similar  results  were 
obtained  with  the  Moody  potential  for  Ni  as  well  as  with 
the  Mishin  potential  for  Cu.  In  these  simulations,  the  dis¬ 
sociation  of  the  screw  dislocation  either  on  the  (111)  glide 
plane  or  on  the  (111)  cross-slip  plane  results  in  extra  ledge 
surface  area  formation  at  the  free  surfaces.  However,  the 
fully  glide  plane  as  well  as  the  fully  cross-slip  plane  core 
structures  are  completely  constricted  at  one  of  the  free  sur¬ 
faces  (from  Escaig  stresses  induced  by  image  stresses), 
whereas  the  partially  cross-slipped  core  structures  are 
spread  out  at  both  free  surfaces.  Therefore,  the  energy  of 
the  partially  cross-slipped  core  structures  relative  to  the 
fully  glide  plane  or  fully  cross-slip  plane  core  structures, 
dy,  can  be  approximately  written  as  dy  =  0.5{yidb  E^), 
where  is  the  formation  energy  of  a  negative  constriction 
[2],  b  is  the  Burgers  vector  of  the  screw  dislocation,  d  is  the 
Shockley  partial  spacing  of  the  partially  cross-slipped  core 
structure  at  the  free  surface  and  71  is  the  free  surface  ledge 
energy  per  unit  area.  Taking  d^^^A  —  5b  [7],  and  the  nega¬ 
tive  constriction  formation  energy  as  1.23eV  [7]  for 
the  Mishin  potential  for  Ni,  and  with  dy  determined  as 
-0.06  eV  from  the  simulations,  one  obtains  the  free  surface 


Fig.  1.  A  3-D  dislocation  line  representation  of  the  fully  glide  plane  structure  formed  for  a  y2[l  1 0]  screw  dislocation  intersecting  two  free  surfaces  at  right 
angles.  The  screw  dislocation  was  formed  with  a  (1 1 1)  glide  plane  cut.  The  x-axis  is  along  [1-10],  the  y-axis  is  along  [1 1 2]  and  the  z-axis  is  along  [1 1 1].  The 
stacking  fault  region  is  shaded  in  pink.  The  colored  arrows  denote  the  Burger’s  vector  magnitude  and  direction  of  the  Shockley  partials.  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  2.  A  3-D  dislocation  line  representation  of  the  partially  cross-slipped  structure  formed  for  a  V2[l  1 0]  screw  dislocation  intersecting  two  free  surfaces  at 
right  angles.  The  screw  dislocation  was  formed  with  a  (1 1 1)  glide  plane  cut.  The  x-axis  is  along  [1 1 0],  y  axis  is  along  [1 1 2]  and  the  z-axis  is  along  [111]. 
The  stacking  fault  region  is  shaded  in  pink.  The  colored  arrows  denote  the  Burger’s  vector  magnitude  and  direction  of  the  Shockley  partials.  (For 
interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


ledge  energy  71  as  550-675  mJ  m~^.  This  value  is  approxi¬ 
mately  one-third  the  (111)  surface  energy  given  by  the 
Mishin  potential  for  Ni  [12].  Using  the  Mishin  potential 
for  Cu,  —1.3  eV,  dy  ^  0.01  eV,  6b  [6,7]  gives  the 
free  surface  ledge  energy  71  as  538  mJ  m“^,  which  as  before 
is  approximately  one-third  the  (111)  surface  energy  given 
by  the  Mishin  potential  for  Cu  [13].  With  the  Moody 
potential  for  Ni,  ~  — 1.58  eV,  dy  ^  0.14  eV,  6b  [5- 
7]  gives  the  free  surface  ledge  energy  71  as  555  mJ  m“^, 
which  is  also  approximately  one-third  the  (111)  surface 
energy  given  by  Ni  Moody  potential  for  Ni  [11].  We  will 
discuss  below  the  application  of  the  nudged  elastic  band 
(NEB)  method  [7,18]  with  the  initial  and  final  states  as 
the  fully  glide  plane  and  partially  cross-slipped  structures 
of  Fig.  3a-d  to  obtain  the  activation  energy  for  surface 
cross-slip  nucleation  in  fee  Ni  in  the  absence  of  ledge  anni¬ 
hilation  forces. 

3.1.  Screw  dislocation  intersecting  the  free  surfaces  in  an 
inclined  fashion 

In  the  simulation  results  to  be  described  in  this  section 
the  screw  dislocation  was  initially  inserted  in  the  simulation 
cell  using  its  anisotropic  elasticity  displacement  field  with  a 
(111)  glide  plane  cut.  For  the  z-axis  oriented  along  one  of 
the  multi-slip  or  single  slip  orientations,  [100]  or  [011]  or 
[111  ]  or  [4 1 3],  the  screw  dislocation  intersects  the  free  sur¬ 
faces  in  an  inclined  fashion.  Figs.  4  and  5  depict  two  differ¬ 
ent  core  structures  obtained  for  the  screw  dislocation 
intersecting  (010)  free  surfaces  at  45°  using  the  Mishin 
potential  for  Ni  both  in  the  [111]  and  [1 1  —2]  projections 
as  well  as  the  3-D  dislocation  line  representation.  The  sim¬ 
ulation  cell  used  to  obtain  the  core  structures  shown  in 
Figs.  4  and  5  had  the  z-axis  oriented  along  [1 00],  the  x-axis 
along  [0 1 0]  and  the  y-axis  along  [001]  with  fixed  boundary 
conditions  along  “y”  and  “z”  and  free  surface  boundary 
conditions  along  “x”.  Figs.  4  and  5  clearly  show  the  rota¬ 
tion  of  the  screw  dislocation  to  minimize  its  energy  due 
to  Fothe’s  forces  [15-17].  The  screw  dislocation,  as  it 


rotates,  decreases  its  line  length,  whereas  its  line  energy 
per  unit  length  increases  since  it  moves  away  from  the 
screw  orientation.  This  results  in  an  equilibrium  line  orien¬ 
tation  for  the  screw  dislocation  where  its  total  line  energy  is 
a  minimum.  Note  that  for  this  orientation  of  the  simulation 
cell,  the  expected  rotations  on  both  the  (111)  glide  plane 
and  the  (1 1  T)  cross-slip  plane  are  identical  at  equilibrium. 
For  one  core  structure,  the  screw  dislocation  completely 
resides  on  the  (111)  glide  plane  whereas  for  the  other  core 
structure,  the  screw  dislocation  is  partially  cross-slipped 
onto  the  (111)  cross-slip  plane  with  the  formation  of  a  neg¬ 
ative  constriction.  For  both  core  structures,  there  is  ledge 
annihilation  at  the  negative  “x”  surface  and  ledge  forma¬ 
tion  at  the  positive  “x”  surface  due  to  the  rotation.  The 
energy  of  the  partially  cross-slipped  core  structure  is  lower 
relative  to  that  of  the  fully  glide  plane  core  structure  by 
~1  eV.  If  the  Burgers  vector  of  the  screw  dislocation  is 
reversed,  then  the  formation  of  a  partially  cross-slipped 
core  structure  requires  the  screw  dislocation  to  be  on  the 
cross-slip  plane  on  the  negative  “x”  side  and  on  the  glide 
plane  on  the  positive  “x”  side,  which  requires  extra  ledge 
formation  on  both  surfaces  and  is  energetically  prohibitive. 
Therefore,  depending  upon  the  sense  of  rotation,  it  is  ener¬ 
getically  favorable  to  form  the  partially  cross-slipped  core 
structure  for  only  one  sense  of  the  screw  dislocation  Bur¬ 
ger’s  vector,  either  positive  or  negative.  Similarly,  the  fully 
glide  plane  as  well  as  the  partially  cross-slipped  core  struc¬ 
tures  having  a  negative  constriction  were  observed  for  sim¬ 
ulation  cells  with  the  z-axis  oriented  along  [4 1  3],  [0 1 1]  and 
[111].  In  these  simulations,  the  free  “x”  surfaces  were  cho¬ 
sen  such  that  the  amount  of  rotation  of  the  screw  disloca¬ 
tion  due  to  Fothe’s  forces  were  identical  on  the  (1 1 1)  glide 
plane  and  the  (111)  cross-slip  plane.  In  all  cases,  the  par¬ 
tially  cross-slipped  core  structure  formed  for  only  one  sense 
of  the  screw  dislocation  Burgers  vector  and  had  a  lower 
energy  as  compared  to  the  fully  glide  plane  core  structure. 
The  fact  that  there  is  a  local  minimum  in  energy  at  both  the 
fully  glide  plane  state  as  well  as  the  partially  cross-slipped 
state  suggests  that  there  is  a  finite  activation  barrier  for 
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(a) 


(b) 


(c) 


Fig.  3.  A  3-D  dislocation  line  representation  of  the  fully  glide  plane  structure  (a),  fully  cross-slip  plane  structure  (b),  partially  cross-slipped  structures  1 
and  2  (c  and  d)  formed  for  a  V2[l  1 0]  screw  dislocation  intersecting  two  free  surfaces  at  right  angles.  The  screw  dislocation  was  formed  with  a  (00 1)  cube 
plane  cut.  The  x-axis  is  along  [1 1 0],  the  y-axis  is  along  [1 1 2]  and  the  z-axis  is  along  [111].  The  stacking  fault  region  is  shaded  in  pink.  The  colored  arrows 
denote  the  Burger’s  vector  magnitude  and  direction  of  the  Shockley  partials.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader 
is  referred  to  the  web  version  of  this  article.) 
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surface  cross-slip  [7],  though  it  is  expected  to  be  small  as 
shown  later.  However,  the  amount  of  ledge  annihilation 
and  screw  dislocation  rotation  may  not  be  identical 
between  the  fully  glide  plane  and  the  partially  cross-slipped 
structures,  and  therefore  the  determination  of  activation 
energy  for  surface  cross-slip  nucleation  becomes  problem¬ 
atic  for  this  geometry. 

4.  Surface  cross-slip  nucleation  activation  energy  in  the 
absence  of  ledge  annihilation  forces 

The  NEB  method  [7,18]  was  used  to  determine  the  acti¬ 
vation  energy  for  surface  cross-slip  nucleation  in  the 
absence  of  ledge  annihilation  forces.  The  fully  glide  plane 
state  and  the  partially  cross- slipped  state  of  Fig.  3a-d  were 
taken  as  the  initial  and  final  states  for  the  NEB  calcula¬ 
tions.  A  total  of  16  replicas  (including  the  initial  and  final 
states)  were  used  in  the  NEB  calculations,  together  with 
an  inter-replica  spring  constant  of  0.04  eV  A~^  [7].  Con¬ 
strained  energy  minimization  for  these  calculations  were 
performed  using  a  damped  dynamics  technique  for  up  to 


20,000  fs  [5].  The  NEB  calculations  were  used  to  determine 
the  energy  profile  as  the  screw  dislocation  transformed 
from  the  fully  glide  plane  state  to  the  partially  cross-slipped 
state.  The  energy  profile  shows  a  maxima  in  between  the 
initial  and  final  states  corresponding  to  the  activated  state. 
The  energy  difference  between  the  activated  state  and  the 
fully  glide  plane  state  (initial  state)  was  taken  to  be  the  acti¬ 
vation  energy  for  surface  cross-slip  nucleation  in  the 
absence  of  ledge  annihilation  forces.  Fig.  6  shows  the 
[111]  and  [1 1 2]  projections  of  the  activated  state  obtained 
using  the  Ni  Mishin  potential.  The  activated  state  has  a 
very  small  portion  of  the  screw  dislocation  on  the  cross-slip 
plane.  The  activation  volume  for  surface  cross-slip  nucle¬ 
ation  with  respect  to  Escaig  stresses  can  be  estimated  from 
the  area  of  this  small  cross-slipped  portion  to  be  extremely 
small.  The  activation  energy  for  surface  cross-slip  nucle¬ 
ation  was  found  to  be  ~0.05  eV  with  the  Mishin  potential 
for  Ni.  Similar  NEB  calculations  were  performed  with  the 
Moody  potential  for  Ni  as  well  as  the  Mishin  potential  for 
Cu  and  the  activation  energy  for  surface  cross-slip  nucle¬ 
ation  was  determined  to  be  ~0.20  and  ~0.09  eV,  respec- 


(a) 


(b)  (C) 


Fig.  4.  A  3-D  dislocation  line  representation  (a)  of  the  fully  glide  plane  structure  formed  for  a  VilllO]  screw  dislocation  intersecting  two  (010)  free 
surfaces  in  an  inclined  fashion.  The  screw  dislocation  was  formed  with  a  (1 1 1)  glide  plane  cut.  The  x-axis  is  along  [0 1 0],  the  y-axis  is  along  [001]  and  the 
z-axis  is  along  [100].  The  stacking  fault  region  is  shaded  in  pink.  The  colored  arrows  denote  the  Burger’s  vector  magnitude  and  direction  of  the  Shockley 
partials.  Also  shown  are  the  [111]  (b)  and  [1 1 2]  (c)  2-D  projections  of  the  position  of  high-energy  atoms  at  the  screw  dislocation  core.  Note  the  rotation  of 
the  screw  dislocation  due  to  Lothe’s  forces. 
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Fig.  5.  A  3-D  dislocation  line  representation  (a)  of  the  partially  cross-slipped  structure  formed  for  a  Vill  1 0]  screw  dislocation  intersecting  two  (0 1 0)  free 
surfaces  in  an  inclined  fashion.  The  screw  dislocation  was  formed  with  a  (1 1 1)  glide  plane  cut.  The  x-axis  is  along  [0 1 0],  the  y-axis  is  along  [001]  and  the 
z-axis  is  along  [100].  The  stacking  fault  region  is  shaded  in  pink.  The  colored  arrows  denote  the  Burger’s  vector  magnitude  and  direction  of  the  Shockley 
partials.  Also  shown  are  the  [111]  (b)  and  [1 1 2]  (c)  2-D  projections  of  the  position  of  high-energy  atoms  at  the  screw  dislocation  core.  Note  the  rotation  of 
the  screw  dislocation  due  to  Lothe’s  forces.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of 
this  article.) 


Fig.  6.  The  [111]  (a)  and  [1 1 2]  (b)  2-D  projections  of  the  position  of  high-energy  atoms  at  a  Vill  1 0]  screw  dislocation  core  in  the  activated  structure.  The 
screw  dislocation  was  formed  with  a  (00 1)  cube  plane  cut. 
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lively,  with  these  potentials.  Note  that  these  are  results  for 
the  activation  energy  for  surface  cross-slip  nucleation  in  the 
absence  of  ledge  annihilation  forces. 

5.  Discussion 

The  results  of  this  work  clearly  show  that  under  certain 
favorable  conditions  (favorable  sense  of  screw  dislocation, 
Lothe’s  rotation  forces  and  minimal  ledge  annihilation 
forces)  surface  cross-slip  can  nucleate  in  fee  Ni  and  Cu. 
This  is  an  important  result  since  it  qualitatively  explains 
the  profusion  of  cross-slip  and  dislocation  multiplication 
in  micron-sized  fee  crystals.  This  finding,  along  with  previ¬ 
ous  findings  on  cross-slip  nucleation  at  screw  dislocation 
intersections  in  fee  materials,  should  allow  higher-level 
mesoscale  models  of  dislocation  behavior  to  better  repre¬ 
sent  the  cross-slip  process  without  resorting  to  ad  hoc 
assumptions  about  obstacles. 

The  activation  energy  for  surface  cross-slip  nucleation  in 
the  absence  of  ledge  annihilation  forces  has  been  deter¬ 
mined  to  be  extremely  low,  0.05  eV  in  fee  Ni  (Mishin 
potential)  and  0.09  eV  in  fee  Cu,  significantly  lower  (by  a 
factor  of  20-30)  than  cross-slip  nucleation  in  isolation  in 
bulk  via  the  Escaig  mechanism  [2,4-7].  Also,  the  surface 
cross-slip  nucleation  energy  is  significantly  lower  than  acti¬ 
vation  energy  for  cross-slip  nucleation  at  mildly  attractive 
screw  dislocation  intersections  in  fee  Ni  and  Cu  [5-7]. 
However,  cross-slip  nucleation  at  mildly  repulsive  screw 
dislocation  intersections  has  been  found  to  be  athermal 
with  zero  cross-slip  activation  energy  [8].  The  present  sur¬ 
face  cross-slip  nucleation  results  are  expected  to  be  useful 
in  understanding  the  mechanical  behavior  of  micron-sized 
pillars  under  pure  tension  and  compression  [19,20].  How¬ 
ever,  it  is  unclear  what  effect  the  presence  of  an  oxide  film 
at  the  surface  has  on  the  cross-slip  nucleation  activation 
energy. 

The  low  activation  energy  for  surface  cross-slip  nucle¬ 
ation  in  fee  Ni  and  Cu,  determined  in  this  work  are  in 
agreement  with  X-ray  topographic  observations  of  prefer¬ 
ential  cross-slip  at  the  surface  in  Al,  Cu  and  Si  [15- 
17,21-23]  and  in  situ  transmission  electron  microscopy 
observations  of  spontaneous  cross-slip  in  Au  foils  [15- 
17,24].  The  low  activation  energy  for  surface  cross-slip 
nucleation  in  fee  Cu  shown  in  the  present  work  is  also  in 
agreement  with  surface  observations  of  cross-slip  of  dislo¬ 
cations  in  alloys  with  very  low  stacking  fault  energies 
(Cu-12-14  at.%  Al)  when  deformed  in  tension  [15-17,25] 
and  when  fatigued  [15-17,26]  despite  the  wide  dissociation 
of  dislocations  (^lOnrn)  and  the  low  temperature  of 
observation  (300  K).  In  addition,  easy  surface  cross-slip 
nucleation  has  been  used  to  explain  the  observed  operation 
of  slip  systems  with  low  Schmid  factors  in  fee  metals  [15- 
17,27]. 

In  this  work,  the  surface  cross-slip  nucleation  activation 
energy  was  determined  in  the  absence  of  ledge  annihilation 
and  Lothe’s  rotation  forces  using  molecular  statics  calcula¬ 
tions.  It  is  suggested  that  if  ledge  annihilation  and  Lothe’s 


rotation  forces  are  present,  surface  cross-slip  nucleation 
will  occur  after  the  screw  dislocation  rotates  and  bends  to 
annihilate  ledges  and  decrease  its  line  energy.  However, 
such  calculations  are  not  amenable  to  molecular  statics 
simulations.  Large-scale  molecular  dynamics  simulations 
at  temperature  may  be  required  to  interrogate  surface 
cross-slip  nucleation  in  the  presence  of  these  forces.  Lor 
some  geometries,  Lothe’s  rotation  forces  may  be  asymmet¬ 
ric  between  the  glide  and  cross-slip  planes  [15-17],  and  in 
such  cases,  it  is  thought  that  the  screw  dislocation  would 
prefer  to  lie  on  a  plane  where  the  decrease  in  line  energy 
due  to  Lothe’s  rotation  forces  is  the  highest. 

6.  Summary 

The  simulations  in  this  work  show  the  nucleation  of 
cross-slip  at  a  free  surface  in  fee  Ni  and  Cu.  The  results 
can  be  summarized  as  follows: 

(1)  Two  different  cores  structures,  one  fully  on  the  glide 
plane  and  the  second  partially  cross-slipped  with 
the  formation  of  a  negative  constriction,  are  found 
for  a  screw  dislocation  intersecting  free  surfaces  at 
right  angles  and  formed  using  a  (111)  glide  plane 
cut.  The  screw  dislocation  is  bent  inwards  at  both  free 
surfaces  to  decrease  its  ledge  energy  in  the  fully  glide 
plane  structure. 

(2)  Lour  different  core  structures — one  fully  on  the  glide 
plane,  one  fully  on  the  cross-slip  plane  and  the 
last  two  crystallographically  equivalent,  partially 
cross-slipped  structures  forming  a  negative  constric¬ 
tion — are  found  for  a  screw  dislocation  intersecting 
free  surfaces  at  right  angles  and  formed  using  a 
(001)  cube  plane  cut. 

(3)  Two  different  core  structures,  one  fully  on  the  glide 
plane  and  the  other  one  partially  cross-slipped  form¬ 
ing  a  negative  constriction,  are  found  for  a  V2[l  —10] 
screw  dislocation  intersecting  two  (0 1 0)  free  surfaces 
in  an  inclined  fashion.  The  screw  dislocation  rotates 
to  minimize  its  line  energy.  The  sense  of  screw  dislo¬ 
cation  has  to  be  favorable  for  the  formation  of  the 
partially  cross-slipped  structure. 

(4)  Surface  cross-slip  nucleation  energy  has  been  deter¬ 
mined  to  be  0.05  eV  in  fee  Ni  and  0.09  eV  in  fee 
Cu,  in  the  absence  of  ledge  annihilation  forces. 
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